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We report a comprehensive study of the paradigmatic quasi-lD compound (TaSe^I performed by 
means of angle-resolved photoemission spectroscopy (ARPES) and first-principles electronic struc- 
ture calculations. We find it to be a zero-gap semiconductor in the non-distorted structure, with 
non-negligible interchain coupling. Theory and experiment support a Peierls-like scenario for the 
CDW formation below Tcdw = 263 K, where the incommensurability is a direct consequence of the 
finite interchain coupling. The formation of small polarons, strongly suggested by the ARPES data, 
explains the puzzling semiconductor-to-semiconductor transition observed in transport at Tcdw- 

PACS numbers: 71.45.Lr, 71.10.Pm, 74.25.Jb, 71.15.Mb 



One-dimensional systems exhibit a rich phenomeno- 
logy due to their reduced phase space. Typically, the 
metallic state is unstable against competing broken- 
symmetry phases such as charge- or spin-density waves 
[TH3]. In the dominant Peierls scenario, these transitions 
are described as electronic instabilities driven by the nes- 
ting properties of the Fermi surface. Alternative models 
have also been proposed [4-6 , and recently even the re- 
levance of the Peierls scenario for real materials has been 
questioned [3 [8]. 

The chain compound (TaSe^I is a paradigmatic 
quasi-one-dimensional (ID) material [9 . At Tcdw = 
263 K it undergoes a charge density wave (CDW) transi- 
tion, accompanied by an incommensurate structural dis- 
torsion [TOj [TT]. As expected, the instability affects the 
electronic structure, and clear signatures of the transi- 
tion are found in the electrical resistivity [I2j fl3] . the 
magnetic susceptibility [14] and the optical response [15] . 
The CDW phase is rather unusual. The main part of the 
lattice modulation is acoustic and transverse [TT], with 
only a much smaller optical component along the chains 
[T6] corresponding to Ta-tetramerization modes [17] . Be- 
sides, structural investigations could not establish the 
existence of a soft phonon mode [18]. The electrical re- 
sistivity is also puzzling. At low temperature, it exhibits 
an activation gap Alt ~ 0.3 eV [13 , but the tempe- 
rature dependence is that of a semiconductor also above 
Tcdw, albeit with a smaller Aht ~ 0.2 eV [19]. The ratio 
^lt/^b^cdw ~ 12, much larger than the 3.5 mean- field 
value, can be rationalized as an extreme manifestation of 
ID fluctuations [20] or, alternatively, of strong electron- 
phonon coupling [18] . 

The electronic structure of (TaSe^I has been explo- 
red both by angle-integrated and by angle-resolved pho- 
toemission spectroscopy (ARPES). Early results repor- 
ted the absence of a metallic Fermi edge in spectra of 



the high-temperature phase [2TJ [22] , suggesting the pos- 
sible manifestation of the Luttinger liquid behavior [23J. 
It was later recognized that strong electron-phonon cou- 
pling leads to the formation of polaronic quasiparticles, 
and to a dramatic suppression of the coherent spectral 
weight [24]. Further ARPES data revealed signatures 
of inter-chain coupling [25] and of the incommensurate 
CDW periodicity [26]. 

In this Letter, we present an extensive survey of k- 
space by ARPES, and first-principles calculations of the 
band structure of (TaSe^L The combination of theory 
and experiment provides an unprecedented view of the 
three-dimensional electronic structure of this compound. 
Surprisingly, (TaSe^I is found not to be a metal, but 
a zero-gap semiconductor. The incommensurate perio- 
dicity of the CDW appears as a direct consequence of 
the finite interchain coupling. The ARPES spectral line 
shapes, typical of small polarons, suggest a natural expla- 
nation for the puzzling semiconductor-to-semiconductor 
transition observed at Tcdw in transport data. 

ARPES experiments were performed at beamline 7.0.1 
of the Advanced Light Source (Berkeley). The energy and 
angular resolution of the hemispherical Scienta R4000 
analyser were set to 25 meV and 0.1°, respectively. The 
polarization and scattering geometry is described el- 
sewhere [27]. We used (TaSe^I single crystals grown 
by the chemical transport method, with typical sizes 
5 x 0.2 x 0.2 mm 3 . They were cleaved in situ at 100 K 
and a pressure of 10 ~ 10 mbar, and measured at 100 K, in 
the CDW phase. First-principles electronic structure cal- 
culations were performed within the density functional 
theory (DFT) framework using the generalized gradient 
approximation (GGA). Spin-orbit effects were included 
by means of the fully relativistic norm-conserving pseu- 
dopotentials acting on valence electron wave functions 
represented in the two-component spinor form [28 . We 
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Figure 1: (color online), (a) Crystal structure of no n- distorted (TaSe^I showing the natural cleavage plane (pink). The 
conventional tetragonal unit cell is outlined, (b) A projection of the structure onto the cleavage plane. t\\ and t± are the 
intrachain and interchain nearest-neighbour hopping integrals in the tight-binding model discussed in the text. ARPES intensity 
maps illustrate the band dispersion (c) parallel (along XP) and (d) perpendicular to the chains, for k\\ — 0.75 tt/c. (e) 3D 
BZ. (f),(g) Constant energy maps extracted at the top of the Ta 5d z 2 band (Eb = 0.3 eV, T = 100 K) in two perpendicular 
planes, as depicted in the 3D BZ. Photon energy in (f) is set to hv — 88 eV, that corresponds to the dashed line in (g). Second 
derivative is used in (g). 



used the Quantum-ESPRESSO software package [29] . 

Figure [TJa) shows the crystal structure of (TaSe^I 
with its conventional tetragonal unit cell (a = b = 
9.531 A and c = 12.824 A). The Ta atoms form chains, 
surrounded by Se4 rectangular units and separated by 
I - ions. Two adjacent Se4 rectangles are rotated by 45°, 
leading to a periodicity c = 4<ixa-Ta along the chain. 
The Bravais lattice actually is not tetragonal, but body- 
centered tetragonal (bet, space group 7422 [9]). The cor- 
responding rhomboedral unit cell contains two formula 
units (4 Ta atoms), instead of four (8 Ta atoms) for 
the conventional cell. The iodine ions define the natu- 
ral cleavage plane, perpendicular to e n = (a + b)/y/2a. 
The in-plane unit vectors parallel and perpendicular to 
the chains are ey = c/c and e± = (a — b)/y/2a, res- 
pectively. The (ej|,e_j_) and (en,e n ) planes are equivalent, 
reflecting the four- fold symmetry of the space group. Fi- 
gure [TJb) shows the projection of the structure onto the 
(ey ,e_j_) cleavage plane. The interchain distance is a/y/2 = 
6.739 A. Figure [jje) shows the bet BZ, which exhibits 
four- fold symmetry, with equivalent k n and k± directions. 
Reciprocal distances are TZ = 2tt/c = 0.490 A -1 (ins- 
tead of tt/c for an isolated chain), XP = tt/c = 0.245 A -1 
and TX = ir/(a/V2) = 0.466 A -1 . The CDW exhibits 
eight equivalent domains with incommensurate wave vec- 



tors 5q = (±a, d=a, ±7), where a — 0.045(27r/a) and 
7 = 0.085(2tt/c) [18]. 

The in-plane band dispersion measured by ARPES 
along and perpendicular to the chains is shown in 
Figs.[TJc) and[TJd). Along the chains, the XP direction 
is probed at this photon energy (hv = 88 eV; see below). 
We observe the characteristic "V-shaped valence band 
(VB) [26 with a minimum at ku = 0. The band maxima 
are at kn — ±tt/c = 0.245 A -1 and 0.3 eV binding energy. 
In a simple picture, it represents the fully occupied, lo- 
west of four subbands split from the nominally quarter- 
filled band, after gaps are opened at the BZ boundary 
by the x4 in-chain periodicity [30]. This band is deri- 
ved from Ta bd z 2 orbitals which overlap strongly along 
the chain, yielding a small effective mass (0.3 m e , where 
m e is the bare electron mass) and an overall band width 
W\\ ~ 4 eV. Strong polaronic effects drastically reduce 
the coherent quasiparticle (QP) weight. The largest in- 
tensity in the ARPES map does not coincide with the QP 
energy, but rather with the maximum of the incoherent 
part of the spectral function. Its energy is lower than the 
QP energy by the polaron binding energy E p w 0.2 eV 
[24] . Weaker replicas of the VB are visible at larger wave 
vectors, in higher-order BZs. 

The VB also disperses in the perpendicular (fej_) di- 
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Figure 2: (color online). Schematic constant energy surfaces 
from (a) a 2D and (b) a 3D tight-binding (TB) model, (c-d) 
Comparison of experimental CE contours measured at Eb — 
0.6 eV (hv — 88 eV), and theoretical contours derived from 
a 3D TB model. Projected 3D (solid line) and 2D (dashed) 
BZs are also indicated 

rection, with periodicity 2n/(a\^2) (Fig.JIJd)). The band 
width is W± « 0.2 eV, corresponding to a smaller but 
finite inter-chain coupling. This is further illustrated by 
a constant energy (CE) map, measured at the top of the 
band (Fig. [jjf)). The open contours are the hallmark 
of a ID system, while the periodic ondulations reveal, 
once again, the interchain coupling [25]. The overall fea- 
tures of the experimental dispersion can be reproduced 
by a simple nearest-neighbour tight-binding (TB) model. 
The hopping parameters can be estimated from the band 
widths. On a rectangular lattice, Wt u±) = 4£y(j_). We ob- 
tain t± « 0.05 eV and t\\ « 1 eV, and (t\\/t ± ) 2 « 400 in 
remarkable agreement with the anisotropy of the electri- 
cal conductivity [3TJ [32] . 

Figure [TJg) reveals a previously unreported dispersion 
along the k n direction, perpendicular to the surface. The 
CE map was extracted from a data set where the pho- 
ton energy was varied between 80 and 150 eV. The inner 
potential used to determine k n was Vq = 16 eV. We ob- 
serve the same wiggling contours as in the (fey, k±) plane, 
consistent with the four-fold symmetry of the lattice. As 
anticipated, the line in k space probed at hv = 88 eV 
(white dashed line) almost coincides with the XP direc- 
tion in the BZ. The separation between the two contours 
is minimum ("waist") along this line, and maximum 
("belly"), along the YZ direction. This difference in the 
electronic structure between YZ and XP plays an im- 
portant role in the CDW formation, as shown below. 

In ID CDW materials, it is often necessary to consi- 
der transverse coupling to reproduce the periodicity of 



the distorted phase. Typically, only one perpendicular 
direction is relevant. Such an effective 2D scenario, sche- 
matically illustrated in Fig. [2|a) by CE surfaces obtai- 
ned from a 2D TB model, is inadequate for (TaSe^I. 
The two inner sheets correspond to the main band, while 
the outer two are folded replicas. The (fc n ,fcj_) planes 
at fcy = ±7r/c define mirror planes between the main 
band and the replicas. Furthermore, (fcy , k±) CE contours 
(black lines) do not depend on k n . By contrast, the expe- 
rimental CE contours of Fig. [2|c) exhibit mirror planes 
at ±27r/c = ±0.490 A -1 . An agreement with experiment 
is found for a 3D TB model [Fig. ^b)]. The constant 
energy sheets are now warped in both the e± and e n 
directions, and the CE contours vary with k n . Indeed 
from k n = (red) to 7r/(a/\/2) (green), even if the 
shape of the contours does not change, the kn separa- 
tion between opposite branches increases. The calcula- 
ted contours for k n = 7r/(a/\/2) (corresponding to XP, 
Figure |2jd)) clearly exhibit the periodicity of the projec- 
ted 3D BZ (solid line) and not that of the 2D BZ (dashed 
line) . Consideration of the full 3D band structure is there- 
fore essential for a proper description of (TaSe^I. This 
observation provides a rationale for the low anisotropy 
of the critical fluctuations measured by neutrons above 
^cdw PS] • 

First-principles calculations, performed for the undis- 
torted (high-T) structure, provide important clues on the 
nature of the instability. The band dispersion is shown in 
Fig. |3ja) for the main high-symmetry directions. Along 
the XP direction, parallel to the chains, two d z 2-\ike 
bands (VB and CB), with extrema at P (fcy = 7r/c), 
are separated by a small gap of 50 meV. The Fermi 
level intersects VB slightly below the top of the band, 
yielding a narrow hole pocket. By contrast, along the pa- 
rallel YZ direction the gap between VB and CB collapses 
at Ep. The top of VB and the bottom of CB merge at 
Ep, yielding a point-like Fermi surface, reminescent of 
the Fermi surface of graphene. Therefore, the band struc- 
ture along YZ is that of a zero-gap semiconductor. The 
crossing point is not at n/c, but at the incommensurate 
kp ~ 1-1 tt/c, such that 2kp « q^ DW • This discrepancy, 
which reflects the small but noticeable interchain disper- 
sion, has important consequences on the properties of the 
CDW instability. 

Figure |3ja) illustrates the calculated Fermi surface 
(FS) of (TaSe^I. The Fermi level was shifted just above 
the bottom of the conduction band to account for a slight 
n-doping which is always observed in real samples as a 
result of iodine vacancies. The FS consists of pairs of very 
flat 'pancakes', perpendicular to the chain direction. The 
largest pair, along TZ, is nested by the incommensurate 
CDW wave vector g^ DW , consistent with a Peierls scena- 
rio. To check that hypothesis, we have performed DFT 
calculations for the distorted structure, where Ta atoms 
are mainly displaced along the chains [16] [17] . For prac- 
tical reasons, we considered a commensurate, tetrameri- 
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zed structure (gy = 2tt/c). The computed band disper- 
sions along TZ in the high-T (red) and low-T (green) 
structures are compared in Fig.[3jc). As expected, a gap 
(Acdw = 0.2 eV) opens between the valence and conduc- 
tion bands. Remarkably, the largest energy gain is not at 
7r/c, but at k = 0.5 x q^ DW > tt/c. This observation 
suggests that the largest gap, and therefore the largest 
electronic energy gain, is achieved for the actual incom- 
mensurate CDW wave vector. Since the elastic energy 
cost is fairly insensitive to small wave vector changes, 
these results support the picture of a Fermi surface dri- 
ven instability. 

The consistency of this scenario is confirmed by the 
ARPES data of Fig. [3jc), which shows energy disper- 
sion curves (EDCs) measured along the TZ direction in 
the CDW state. The feature associated with VB indeed 
reaches a minimum binding energy for fcjj (blue curve) 
larger than tt/c (dashed). Moreover, a second peak ap- 
pears at lower binding energy in a small k range around 
fcjj. It is the signature of the occupied CB states. The 
experimental CDW gap, defined as the minimum energy 
separation between the two features, is in good agree- 
ment with the theoretical prediction. As for VB, strong 
polaronic effects are expected to alter the spectral weight 
distribution of the CB. The measured peak overestimates 
the actual binding energy by the polaron binding energy 
E p ~ 0.2 eV [24]. The existence of this additional spectral 
feature was already suggested by previous data [24 , but 
its true nature was not recognized because a much more 
limited region of k space was probed in that experiment. 

Finally, the proposed scenario naturally explains the 
puzzling semiconductor-to-semiconductor transition ob- 
served in transport at Tqdw [ISl [19] . For T > Tcdw? the 
activated behaviour (Aht ~ 0.2 eV) reflects the diffusive 
motion of small polarons, which reduces the conductivity 
by a factor e - E ^l k ^ T ^33 For T < t CD w, the activa- 
tion energy is increased by the opening of the CDW gap. 
Transport gives Alt ~ 0.3 eV, that compares well to 
Aht + Aqdw = 0.35 — 0.4 eV. In this perspective, esti- 
mates of the ratio AltAb^cdw an d the importance of 
ID fluctations, that act on Acdw only, should be recon- 
sidered. 

In summary, an extensive survey of k space by ARPES, 
and new first-principles calculations, show that the ID 
compound (TaSe^I is not a metal, but a zero-gap se- 
miconductor with non- negligible transverse interactions. 
Theory and experiment support a Peierls-like scenario, 
where the incommensurate CDW vector is a direct conse- 
quence of finite interchain coupling. By taking the pola- 
ronic nature of QPs, the model also explains puzzling 
properties of (TaSe^I, such as the semiconducting cha- 
racter of the resistivity on both sides of the transition, 
and the unusually large value of the ratio Alt/^b^cdw- 
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Figure 3: (color online), (a) Electronic band structure of non- 
distorted (TaSe^I calculated from first principles. The inset 
shows the Fermi surface. The blue arrow depicts the CDW 
vector, (b) Distorted (green) vs. non-distorted (red) band 
structures along TZ. The horizontal lines indicate the Fermi 
level in the intrinsic (dotted) and n-doped (plain) cases, (c) 
EDCs along the TZ direction (hv — 110 eV, T < T CD w). 
Natural iodine vacancies yield an effective n-doping that al- 
lows to probe CB. Small polaron formation suppresses the QP 
weight at Ef and transfers the spectral weight at a higher E p 
binding energy. 
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